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1 Introduction 



Local boundary conditions of chiral bag type for Euclidean Dirac operators have at- 
tracted, over many years, the interest of both mathematicians and physicists. 

From the physical point of view, such conditions are closely related to those appearing 
in the effective models of quark confinement known as chiral bag models They 
contain a real parameter 9, which is to be interpreted as the analytic continuation of the 
well known ^-parameter in gauge theories. Indeed, for 9^0, the effective actions for 
Dirac fermions contain a CP-breaking term proportional to 9 and proportional to the 
instanton number [21 E|. Moreover, these boundary conditions can be applied to one-loop 
quantum cosmology jU |3],supergravity theories [0] and branes [7J and are important in 
the investigation of conformal anomalies jHj- Applications to finite temperature problems 
were studied in El EH 

From a mathematical point of view, as shown in [3], Euclidean Dirac operators under 
local boundary conditions of chiral type define self-adjoint boundary problems. More 
recently, it was shown that both the first order boundary value problem and its 
associated second order one are strongly elliptic. The asymptotics of the smeared trace 
of the heat kernel was studied, for the second order problem on general Riemannian 
manifolds with boundary, in ^Tj, where use was made of functorial methods [T?? | I14 [ ITK] 
and special case calculations presented in ^Hl HZ! > m this context see also [Tj\ . 

The main characteristic of these boundary conditions, as compared to nonlocal or Atiyah- 
Patodi-Singer JH] ones, is the explicit breaking of chiral symmetry. Since the Dirac 
operator has no zero modes this breaking comes only from the excited modes. Since the 
asymmetry is encoded in the corresponding eta-function, the study of this function if of 
central importance, physical applications for example arising in the analysis of fermion 
number fractionization in different field theory models [T§1 |2*UI |2~H |2~2~] . 

A recent step forward in the analysis of the eta function is PJJ, where the eta function 
associated to the boundary value problem was shown to be regular at s = 0, and 
where some properties of the asymptotic coefficients in the trace of the smeared kernel 
corresponding to the eta function were obtained. 

In the present paper, we evaluate the spectral asymmetry for Dirac operators on the ball 
in two and four Euclidean dimensions, with their domain defined by local boundary con- 
ditions of chiral bag type. Moreover, making use of these results, and the corresponding 
ones for cylindrical manifolds obtained in [T2~l l2~H] we determine some properties of the 
leading asymptotic coefficients in the trace of the smeared kernel corresponding to the 
eta function, for arbitrary Riemannian manifolds. 

The outlay of the paper is as follows: In section |2] we present a new method, based 
on Group Theory, to obtain the spectrum of the boundary value problem under study. 
Section 01 presents an integral formula for the spectral asymmetry of the Euclidean Dirac 
operator in a ball of arbitrary even dimension. Section |U is devoted to the explicit calcu- 



2 



lation of the asymmetry in a disk, while the specral asymmetry in the four dimensional 
case is evaluated in section Finally, in section ||3 some properties of the leading coeffi- 
cients in the heat trace associated to the eta function are established through functorial 
methods, and the use of special cases. To be self-contained, the Appendix contains a 
derivation of the Debye expansion of the Bessel functions. 



2 Eigenvalue Problem for the Dirac Operator 

Although the eigenvalue equations for the problem at hand were derived before |2HEj, 
we present here a different derivation, based on group theory. We write the free Euclidean 
Dirac operator (i.e. gauge fields are absent) as 

where we used the Einstein summation convention for the index A — 1, . . . , D, and the 
Dirac-matrices fulfill the Clifford algebra 

{IAiIb} = 25 AB . (2) 

We have to choose boundary conditions for the spinors such that the Dirac operator 
is Hermitian and in the following we will consider chiral-bag boundary conditions. Given 
the projection operators 

n T =~(lTi7*e> e S) (3) 

with free parameter 9 and 



7* = (-i) d 7i • • • 7£>, S = j A x A /r, r = ^x A x A , (4) 
these boundary conditions are defined as 

n-^U = o. (5) 

In particular, S is the projection of 7 onto the outward unit normal vector. In contrast 
to Atiyah-Patodi-Singer boundary conditions, the chiral-bag boundary conditions are 
local. One can further show that, for simply connected boundaries, these boundary 
conditions do not allow for zero modes. Further details can be found in 

P and Il_ commute with the total angular momentum 

Jab = L AB + Sab, L A b = -i(x A d B - x b 0a), ^ab = ^[jA, 1b\- (6) 
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Using the algebraic approach developed in we can first diagonalize the total angular 
momentum, i.e. determine the spin spherical harmonics by group theoretical methods. 
Our aim is to construct the highest weight states. These states are eigenstates of the 
Cartan operators of the so(.D)-algebra and are annihilated by the corresponding raising 
operators. The remaining states in each multiplet can be obtained by applying low- 
ering operators on these highest weight states. It is appropriate to introduce complex 
coordinates and creation/annihilation operators as follows, 

z a = x 2a -i + ix 2a , d a = ^(<9 a;2a _ 1 - id X2a ), (7) 
ipl = o (72a-i + 172a), *Pa = o(72a-i - 172a), a = 1, . . . , d = D/2. 



For the fermionic operators, one easily verifies the relations 

{ip a , 4} = 5 ab , {^a, M = 0, {tf>l 4} = 0. (8) 

Thus, tp a ijp a ) acts as a creation (annihilation) operator, and we can employ the usual 
(fermionic) Fock space construction, starting with the vaccum state 

V>.|ft> = 0, |ai...a m )=^ 1 ...VLl n ). ±<™<d. (9) 

Now we choose the following Cartan- Weyl basis for the so(D)-algebra: Cartan operators 
H a and raising operators Ei,...,E d (corresponding to simple positive roots) read (no 
sum!) 

H a = Z a 8 a - Z a B a + -(ipl^a ~ i'ai'l), CL=l,...,d, 

E a = -i(z a d a+ i - Za+i&a + iplipa+i), a = 1, . . . , d - 1, (10) 
E d = -\{z d -id d - z d d d „ x + ri-iV'l)- 

Since the operators in (fTUj) act trivially on the radial part of spinor wave functions, we 
consider the angular part only. We can easily determine the highest weight states with 
respect to the orbital part, they are given by 

h = 4 (ii) 

These states are annihilated by all simple positive roots, and the eigenvalues with respect 
to the Cartans read (Hi, . . . , H d ) = (£, 0, . . . , 0). Similarly, for the 'fermionic' part there 
are only two highest weight states given by 

X + =4-'^M and X ~ =4---4-i\ty. (12) 

The corresponding eigenvalues of the Cartan operators read (g, . . . , , |, ±|), respectively. 
Next, we determine the highest weight states of 'fermionic' and 'bosonic' degrees of 
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freedom together. Two highest weight states can be constructed easily; they are just 
given by the tensor products 

4>i = <PiX + and 07 = <f> e x~, (13) 

with eigenvalues of the Cartan operators (£+|,|,...,|,±|), respectively Furthermore, 
we observe that the operator S in 

S = {^ a z a + ^ a z a )/r = S\ S 2 = l, (14) 

commutes with the total angular momentum and therefore maps highest weight states 
into highest weight states. We obtain two additional highest weight states, 

<Fl = S(f)j and 07 = SV>7, (15) 

with eigenvalues §, • • ■ , §, ±|). Equations (fT^j) and (fT5]l contain, in fact, all highest 

weight states, which is seen as follows. The tensor product of the bosonic representation 
(£, 0, . . . , 0) with the fermionic representation (|, . . . , ^) is given by 



(£, 0, . . . , 0) <E> j • • • ) 2) (^~l~2'2' --- '2) ® 2 > 2 > " " " > 2 > 2- 



(16) 

4>t x + <f>t-i 

and the tensor product of the bosonic representation (£, 0,...,0) with the fermionic 
representation (|, . . . , |, — |) is given by 



(£, 0,...,0) (g) (|,...,i, 2) — (^+ 2' I' • • • ' 2' 2) ® 2'2'---'D 

The highest weight states appear at the correct places, which is determined by the degree 
of the polynomials in xa and the chirality of the states. By counting the dimensions of 
the representations and comparing it with Weyl's dimension formula for the groups 
we see that there are no further representations in the tensor product rule. 

Next, we investigate the chiral-bag boundary conditions. As stated above, Il_ commutes 
with Jab, and we can diagonalize Il_ in our basis of highest weigth states. We may 
express 7* with the help of ((Zj) by 

d 

7* = n^»-^)> suchthat 7*l^> = (-) d |fi>. (18) 

a=l 

Since 0/ and 0/ (and likewise 07 and 07) have the same eigenvalues with respect to 
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the Cartan operators H a , we allow for linear combinations of them, 

*t=fi(r)(ft + g?(r)4t. (19) 

Imposing the boundary condition = leads to the following equations for the 

radial functions 

f±(R) T ie ±e gf(R) = 0. (20) 

Finally, using complex coordinates, we solve for the spectrum of the Dirac operator, 

i0Vf = \*f, i$ = 2iij a d a + 2i^ a d a . (21) 

From the eigenvalue equation of the Dirac operator ([21)1. one obtains a system of coupled 
first-order differential equations. Using u — £ + d — 1/2, the system reads 

i(ft)' = \gt, i(gt)' = \ft- 2 jvgi. 
It can be easily solved, and the solutions are given by 

/±(r) = r l ' 2 - v { Cl J^ 1/2 (\X\r) + c 2 N u _ 1/2 (\X\r)), (22) 
gf{r) = -i sign(A)r 1 / 2 -^(c 1 J 1 , +1/2 (|A|r) + c 2 N u+1/2 (\X\r)). (23) 

Finally, keeping only those eigenf unctions that are nonsingular at the origin, imposing 
the boundary condition ()20)1 and defining k = \X\R, we obtain 



Ju-i/2(k) — sign(A)e 61 J v +i/ 2 {k) = for the + case, (24) 
Jv-i/v(k) + sign(A)e~ e J v +x/2(k) = for the — case. (25) 

The degeneracy of the eigenvalues can be determined by Weyl's dimension formula. For 
each £ in D = 2d dimensions and for both, fl23| and (|25|). the degeneracy is given by 
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with d s = 2 d the dimension of the spinor space. 



3 The eta function as an integral for arbitrary D 

Our starting point for the evaluation of the spectral asymmetry will be a contour integral 
representation of the eta function, which we derive first; for the general strategy see, 
e.g., [23 E]- In order to clearly distinguish positive and negative eigenvalues, let us 
write down equations (|2*Ijl and (}25j) for A > and A < 0. 



6 



For A > we have 



while for A < 



J v -i/2{k) - e e J v+1/2 (k) = 0, 
Ju-i/2{k)+e- e J u+l/2 {k) = 0, £ = 0,..,oo, (27) 



Ju-i/2(k) + e e J u+1/2 (k) = 0, 
J„- 1/2 (k)-e- e J u+1/2 (k) = 0, £ = 0,...,oo. (2* 



Hence, the ^-function 

r ? ( S ,i^) = ^(sgnA)|A|- 

A 

is given by the following contour integral in the complex plane, 

( -jA f 1 [a - d i f Ju-i/2(zR)-e e J v+1/2 (zR) \ 

= y,di(D)— dzz -r-logl — -(0^-0), 

^ 2m y dz V Ju-x/2\zR) + e y J u+1/2 (zR) J 

where the contour T encloses the positive real axis counterclockwise. The notation above 
means, that we have to subtract the same expression, with 9 replaced by —9. 

We deform the path of integration such that we integrate along the imaginary axis. After 
using the definition of the modified Bessel functions and the elementary relation [21 

1 1 + xx 

arctan x — — In , 

2i 1 - ix 

we find 

oo 

2 7T s < ~*^ f d 

r](s, 10) = — cos ( — ) > di{D) I dtt~ s — [arctan Q u (9, tR) — arctan Q u {—9, tR)} 
n V 2 / z — 4 I dt 

n n J 



where the notation 



Q v (9,x) = e* L+1/ * {X) 



4-1/2 (x) 

has been used. 

In order to evaluate rj(O,i0), which is by definition the spectral asymmetry, we will use 
in the forthcoming sections the shifted Debye expansion of Bessel functions summarized 
in Appendix El We change the variable in the integral according to t — We observe 
that all terms combine nicely if j2H| 

7r 1 + xy 

arctan x — arctan y — arctan 

2 x — y 



7 



is used. With 

PJ9, x) = 1 ( 4+1 / 2(x) + 7 "-V2(s) 
2sinh0 \I v -i/2(x) I u+1/2 (x) 

we obtain the simple looking result 

rj(s,i$) = — R s cos (— ) y^d e (D)i/- 8 / duu- s — arctan P^w). (29) 
71 V 2 J 7^ J du 



1=0 



Using the notation 



s/1 + u 2 t I 



wsinh# ' y'l + u 2 ' 



as 1/ — > 00, equations (|44|) and ()45|) in Appendix 1X1 show that 

P v (6, vu) = 6(1 + S) 

where 

5 = --t 3 + ^-t 4 (l - t 2 ) + V (l - -i 2 + -t 4 l + (z/- 4 ) . 
Using this in 

b b 3 b 3 3b 2 — 1 

arctan (6[1 + 5]) = arctan b + 5 — 2 - 5 2 — + 5 3 - + 0{5% 

the leading three orders of the u — > 00 expansion are obtained. Introducing Z = 
2 + u 2 + w 2 cosh(2#), where naturally occurs as a factor multiplying 5\ we find 
explicitly 



arctan P u (8, vu) = arctan 



2sinh# 



+ 2u^ r 1 1 , + J_ /5 J3 _ 5 A +1(ti _ 23 (6 + 15 f 
Z [v 2 \8 8 / i> d 

4wsinhfl f J_ 1^ 1 /5^ 4 5^ 



Z 2 [z/ 2 4 z/ 3 

i^r 4 _i Vsinh2 



z/ 3 Z 3 [ 3 

These asymptotic contributions are all we will need for the evaluation of 7/(0, 1$) in two 
and four dimensions. More general, in D dimensions, D — 1 asymptotic orders would be 
needed. 

Note that using (|5Uj) in (J29)l . each order of 1/V leads to the appearance of zeta functions 



8 



of the Barnes type [22 EBl HH , 

(b(s,o) = ( 

In detail we have 



D + £-2 



(£ + a)- 



00 1 / n _ 1 \ 



This fact is characteristic for spectral problems on balls. 

All resulting u-integrations can be done in terms of hypergeometric functions, 

-P 



duu s (l + u 



2\-a 



1 + 



u 



(cosh(2#) - 1) 



2(1 + u 2 ) 

2T(a) 



(31) 



2 Fi ( /3,^-^;a;-sinh 2 # ) . (32) 



This allows us to express all terms resulting from the asymptotic expansions (|50|l in the 
compact form 



A i,k,l,j{s) 



sd. 



cos 



wr(j + ^)r(ti) 



2^'+% v 2 



r (; + t) 



D - 1 



x 



:•/'! ! 7,^;j' + |;-sinh 2 e 



(33) 



We next use these results to evaluate the asymmetry on balls in two and four dimensions. 



4 D = 2: The asymmetry on a disk 

In two dimensions the degeneracy is dg(2) = 1 and we have v = I + 1/2, such that from 
flU we get 

00 

2 7TS f d 

rj(s,i$) = -R s cos( — ) y^~ s / duu s — — arctan P^IO, uv\ (34) 

7r V 2 / ' / cm 

Here, the relevant Barnes boundary zeta function reduces to a Hurwitz zeta function, 
i.e., Cb( s ; \) = Ch{s, 1/2). In this case, in order to obtain the analytic extension to s — 0, 
it is enough to consider only two terms in the Debye expansion. In order to see this we 
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integrate (jHljl by parts and have for < 9ft(s) < 1 



r)(s,i$) 




: s cos[ — ) \ v s / du u s 1 arctanP„(6',Mi/) 



The s-factor in the numerator can be cancelled by singularities coming from divergencies 
in the integral (this is the case for the ^-independent term in the Debye expansion) or 
from the pole in the successive Hurwitz zeta functions (this is the case for the term of 
order v~ l in the Debye expansion). 

The leading two asymptotic terms contributing to 77(0, 1$) are 

T](s, i$) = smh9(A ,i,i,i(s) + Ai, 2 ,i,i(s) + ...) 
Because (h(0, 1/2) = we have Ao,i,i,i(0) = and this shows 



Note that the result is invariant under the transformation 9 — > 6 + lixi. 

The case D = 2 was treated before, in reference as an example of a non-product 
manifold. Unfortunately, the term coming from the pole in the Hurwitz zeta function was 
missing in that calculation, which led to erroneous conclusions about relations between 
results on product and non-product manifolds stated in the same reference. 



In the four-dimensional case we have v = I + 3/2 and the degeneracy reads de(4) = 
v 2 — 1/4. Our immediate concern is the evaluation of the ^-function at the value s = 
in dimension D = 4. Arguing as below ([34)1 . terms up to the order need to be 
considered. Further simplifications occur since (see, e.g., [31H I14j) 




(35) 



5 Asymmetry in D = 4 




As a consequence we find that 



Aw,i(o) 
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for all relevant values of k, Then, the contributions to the eta invariant read 

77(0) = C iAidAk,i,j{0) (36) 

i,k,l,j 

where 2 = 1,3 contributes and the non- vanishing numerical multipliers are found from 
(f3T)|) : in detail, the results are 

23 15 
Ci,2,i,i = sinh6> £3,4,1,1 = 2 sinh 6* £3,6,1,1 — — — sinh^ C 3 ,8,i,i — — sinh 6* 

5 5 1 

£3,4,1,2 = -- sinh 9 6*3,6,1,2 = - sinh 9 £3,4,1,3 = sinh 9 £3,6,3,3 = -- sinh 3 9. 

For all relevant values of i, k, l,j, the hypergeometric functions turn out to be hyperbolic 
functions. Using (see, e.g., (HOI El) 

Res ( B (^s = 3; |^ = ~, Res ( B (^s = 1; = -~, 

and adding up all contributions, we find 
1 tanh ( -) 

77(0, i0) = <^r (259 + 344 cosh 9 + 161 cosh(2#) + 16 cosh(3fl)) . (37) 

6144 cosh (|) 

6 Determination of the leading coefficients in heat 
traces 

In this section, we use the special case of the ball just presented, together with sim- 
ilar results for the cylinder and functorial techniques, to find some results about the 
coefficients a^(F, P, n_) in the expansion [T2"l ITB^ I3T] 



00 

Tr(FPe-* p2 ) ~ ]T <(F, P, n_)r 



n=0 



with P an operator of Dirac type decomposed as 

P = ijjVj + ip. 

Here V denotes a connection on a vector bundle V over a compact D-dimensional 
Riemannian manifold M with smooth boundary dM such that V7 = 0, furthermore 
■0 and F are smooth endomorphisms of V. 

The functorial techniques will consist of conformal transformations and of relations be- 
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tween the above expansion and the well-known expansion for the heat-kernel [T2*| ITH| ITTTj 

oo 

Tr (Fe- tp2 ) ~ ^4(F,P 2 ,n_)^- D ^ 2 . 



n=0 



To write down the geometrical structure of the coefficients a% let L a b be the second 
fundamental form, a,b = 1, ...,P> — 1, and let -d denote the derivative with respect to 
the exterior normal. 

Lemma 6.1 Let f be scalar, and F = f ■ Idy. There exist universal constants di(9,D) 
such that 

a2(F,P,U_) = 0, 

a?(F,P,II_) = (4tt)- d / 2 | (1~D)J fTr(1>)dx + J d 1 (6,D)f Tt{Id)dy 1 , 
a^(F,P,n_) = (4tt)-^- 1 )/ 2 y Ti (d 2 (6,D)L aa f + d 3 (9,D)f^(iS) 

dM 

+d i (9, D)fip(iS) + d 5 (9, D)/VT* + 4(0, + d 7 (e, D)f, D Id v ) dy 

Proof: This follows from the theory of invariants as described e.g. in |13| . □ 

We have used the invariant iS, with S the projection of r y A onto the outward unit normal 
vector, such that the numerical multipliers remain real as is commonly chosen. We next 
determine the universal multipliers di(9, D), i = 1, 7. We first exploit known special 
cases. 

Lemma 6.2 We have 

di(6,D) = i(D - 1) sinh0 2^^1,1-^;?;- sinh 2 
d 2 (6,D) = -l(D-2)sinh^ 2 F 1 ^l,^-^;2;-sinh 2 ^. 

Proof: This follows from the calculations on the ball with / = 1 and if) = 0, which is the 
situation considered in the previous sections. For convenience, we also put the radius of 
the ball R = 1. We will use the relation O ED] 

<(l,P,n_) = ( - + 1 ) Res rj(D-n,P) (38) 



2 V 2 

and, therefore, we can use the previously analyzed //-function r)(s,i$). 

Instead of looking at s — 0, the information relevant for Lemma 16.11 is found by con- 
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sidering r](s, i$) in eq. (|2T)|) about the points s = D — 1 and s = D — 2. Arguing as 
before, we need to consider the leading two asymptotic contributions in ()30|) only and 
find, modulo terms that do not contribute to the relevant residues 

v(s, i$) ^ d ^\tl) sinh 6 2Fl ( 1; I - sinh2 e ) Cb ( s; ^Y^) 

+ ld s s(s + l)sinhfl 2 F 1 ( 1, 2; - sinh 2 9 J ( B + 



The residues of r](s, i$) at s — (D — 1) and s = (D — 2) follow immediately from the 
residues of the Barnes zeta function, see e.g. [32], 

/ D — 1\ 1 1 

Res (b[D-1; 



2 J (D-2)\ r(D-l)' 
= 0. 

Using equation (|38|). and the fact that the volume of the (D — l)-dimensional unit sphere 
is 27r D / 2 /r(D/2), we conclude the proof of Lemma f6. 21 □ 

We next apply conformal variations of Dirac type operators. 

Lemma 6.3 We have 

Proof: Let / be a smooth function with f\dM = 0. Define g^e) := e 2e ^ g^, u (0) and 
P(e) := e~ e fP. Let V be the standard spinor connection. We expand P = vy v Vd v + ip 
with respect to a local coordinate system x = (xi, ...,xd) and use the metric to lower 
indices and define j„. Then the connection transforms like 

V(e)a M :=Va, + ^e(-/ ; , 7 ^ + 

Furthermore, 

m = - ~e(£> - l)fAh u ) 

Note that the boundary condition remains unchanged under conformal variation. Pro- 
ceeding as for the heat kernel coefficients, one shows that the eta invariant coefficients 
satisfy the equation [33] 



d_ 

de 



al (1, P(e), n_) = (D - n)al (/, P, IL_) . (39) 

e=0 
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To study the numerical multiplier dr(9, D) we need the variation of 

L aa = ~fL aa + (D-l)f ;D . 



d 

de 



e=0 

Applying equation shows the assertion. □ 

In order to determine the numerical multipliers d$(9, D), d±(9, D), d§{9, D) and d^(9, D) 
we relate the eta invariant to the zeta invariant. We will then evaluate the zeta invariant 
on the -D- dimensional cylinder for the case of an arbitrary endomorphism valued P. 

The result we are going to use is the following: 

Lemma 6.4 Let F e C°°{End(V)) and let P(e) := P + eP. We then have 
d t al(l, P(e), n_) = (n - D)ai_ 1 (F, P(e), IL_). 



Proof: The proof is insensitive to the boundary conditions imposed and parallels the 
proof in (Ml ESI- □ 

Remark 6.5 The very useful property of this result is that the coefficient for the 
eta invariant is related to the coefficient a^ l _ 1 for the zeta invariant, which will have a 
significantly simpler structure. 

In order to apply Lemma 16.41 to the coefficient a\ we need the general form of the a\ 
coefficient. 

Lemma 6.6 Let F £ C°°(End(V)). There exist universal constants fi(9,D) such that 
ai(F,P 2 ,n_) = 

(47r) -(D-i)/2 J Tr {/s( ^ D ) F% ( iS ) + / 4 (0, D)F(iS) + f 5 (9, D)F % + h{9, D)F} . 

dM 



Proof: This follows immediately from the theory of invariants taking into account that 
F is in general a matrix-valued endomorphism. □ 



Remark 6.7 Lemma \6.4\ relates the universal constant dj(9,D) with fj(9,D), j = 
3, 6. In detail we have 

d j (6,D) = -(D-2)f j (e,D), J = 3,.. .,6. 

Finding the fj(9, D) is easier, because they follow from the case with if) = 0. They can be 
evaluated from the cylinder where the heat kernel is known locally \lty and its coefficients 
can be evaluated for an arbitrary endomorphism P. 

In order to summarize the results of ^2] we need to introduce the relevant notation. Let 
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M = R + x N be an even dimensional cylinder equipped with the metric ds 2 = dx 2 D + ds 2 N , 
where xn is the coordinate in R + and plays the role of the normal coordinate, and 
ds 2 N is the metric of the closed boundary N. The coordinates on N are denoted by 
V — {yij V2, •••) Vd-i)- To write down the heat kernel on M for P 2 = (ijj'Vj) 2 with 
boundary condition n_, we call (j) w (y) the eigenspinors of the operator B = j^SjaVa, 
corresponding to the eigenvalue u, normalized so that 

with S " 1 the Dirac delta function, and 

dy <f>*M<f>M = 1 • 



A' 



Finally we need x = (i/,xd), £ = %d — x' D , i] = xd + x' D , u u (rj,t) = -fj= + Vru;tanh0, 
and the complementary error function 



erfc(x) = —= I d£e ^ . 

'IT 



We then have the result. 

Lemma 6.8 From the calculation on the cylinder M = R + x N we obtain the following 
values for the multipliers fi(9, D), i — 3, 4, 5, 6: 

f 3 (9,D) = ^cosh D - 2 £, 

h(B,D) = 0, 

f 5 (8,D) = - cosh D ~ 2 9 sinh 9, 

f 6 (6,D) = ^(cosh D -^-l). 

Proof: In ^2] we have shown that the local heat-kernel on the cylinder reads 

K(x,x';t) = *Y,<l>M)Mv)e- u * Ue=£ - 1 

+ u2/m 1 ~ V^^tanhge^fa^erfcK^,^) e=^\, (40) 
cosh (9) L J J 

with n + = (1/2) (1 + ie 6 ' 7 *7*S'). (Note that [12] and the present article use different 
conventions. Here we use the exterior normal contrary to the interior normal there. 
Furthermore, the 7* here is minus 7 there. As a result, in the solution formula from 
we have to replace 9 by —9 in order to find a solution for the problem considered 
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here.) The first term is the heat-kernel of the manifold R x N, a manifold that has 
no boundary. Therefore we do not consider this term further as it provides no relevant 
information for Lemma f6. 11 

It is natural to introduce the heat kernel KB(y,y';t) of the operator B 2 , 

furthermore, to make the single steps easier to follow we use the splitting 

V47rr — 



1 



cosh 2 (#) 
7rf)u; tanh Be u - ^ erfc {u„ (rj,t)) 



e 4t . 



We are interested in the trace Tr L 2 (FK(x, x;t)). We assume F = F(y) to be independent 
of the normal variable x&, such that the xu-integration of the L 2 -trace can be done 
without greater complication. 

First, it is straightforward to see that 



dx D F(y)K(x,x;t) = --F{y)K B {y,y;t). 



(41) 



The representation for K 2 (x,x;t) can be conveniently rewritten as to perform the red- 
integration. We have 



K 2 (x,x;t) 



2 cosh 2 9 
d 



e-^ +u * (2w) erfc M2zd, £)) 



dx 



where we used the relation 

1 d 



2dx D 



D 



e~ x y t+u " (2 ^' t) erfc(^(2a; D , t)) 




- u tanh 9 e ul {2x ° ' 4) erfc (u„ (2x D ,t)) 
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Therefore, 

oo 

f dx D FK 2 (x,x-,t) = 2 cq ^ 2 C W^K^n+n; 

e ^ 2 tanh 2 e erfc ^^tanh^ . 

We need to collect the contributions to the coefficient a\(F, P 2 ,n_). The first term, 
equation (j41j) . can be described by the heat-kernel of the boundary and we find the rele- 
vant contribution to be — (l/4)aQ(F, B 2 ). In order to find the contribution of K 2 (x, x; t) 
is considerably harder; we found it most convenient to relate the heat-kernel coefficients 
to the zeta function, 



C(s;F,P 2 ,n_) = f^Jdt t^Trv (Fe~ tp2 ) 

o 



in the standard fashion 



Res C (z; F, P 2 , H_) = fl ^(^ n -) > (42) 



for 2 = P/2, (P - l)/2, 1/2, -(2n + l)/2, n G N. 
For K 2 (x, x; t) the related zeta function contribution is 



oo 

J dt t s ' 1 e"^~e (l + erf(-V*wtanh0)) I . 



o 



The integral can be performed in terms of a hypergeometric function to read 



X 



cosh 2s 



2* 



2 _ / 1\ „ /l 13 . , , > 



r(s) ^=r ( s + - j sinh0 sgn(^) 2 F X I -,s + -; -; - sinh 
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In terms of the boundary spectral functions this is 

cosh 2s ~ 2 9 

C 2 (,;F,P 2 ,n_) = C °^ (g) U {T(s)C(s;U + UlF,B 2 ) (43) 
- A r ( s + sinh 2 F X Q, s + ~; | - sinh 2 ^ ^(2a; n+IT + F, 5) 

In order to find the heat-kernel coefficient af, we use the relation 

/D-i o \ 4(F,p 2 ,n_ 
Res C ;F,P 2 ,n_ ' - 1 v ' ' 



2 ' ' r(^i) ' 

The eta invariant in does not contribute as r/(£) — 1; H+II+F, B) = 0. Therefore, 

fD-l 9 \ 1 n , an (n+n+F, P 2 ) 
Res C 2 (^-g-j F, P 2 , IT_ J = - cosh D " 3 9 ° 1 r +^j >- 

and the contribution to the heat-kernel coefficient a{(F,P 2 , LL) is 

icosh D - 3 fl (47r)-^ ^ Tr(n+n;P). 

N 

To compare it with the form given in Lemma f6. 61 we use 

n + n^ = - cosh 9 (cosh # + 7* sinh # + 7* (iS) ) , 
providing the relevant heat-kernel contribution in the form 

^-cosh D ~ 2 # (47r)-^ y Tr (cosh 0F + sinh 07*F + 7*(LS)F). 

iV 

Adding the contribution from Ki(x,x;t) shows the result. □ 



7 Conclusions 

In this article we have determined the asymmetry r\ (0, i$) of the Dirac operator with 
chiral bag boundary conditions given by (0) on the two-dimensional and four-dimensional 
ball, see equations (J35J) and (pT7j). Furthermore, the leading coefficients in the trace of 
the smeared kernel corresponding to the eta function were obtained, see Lemma 6.1, 6.2, 
6.3 and 6.8. 
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A Debye expansion of Bessel functions 

Consider the ordinary differential equation 

(<9 2 2 + ud z + v)(f)(z) = 0. 
For the solution we use the ansatz 



(f)(z) = exp |y dt p(t) | ib(z). 
With the choice p = —u/2, the resulting differential equation for ip is 

/ 2 

u u 

ib" + qib = 0, where q — v -. 

2 4 

It is more convenient to consider 

d z In (b = p + S, where S — d z In ib. 
Contact between S and the solution is made by observing that 

4>{z) = const exp < I dz p{z) > exp < I dz S{z) 



The differential equation for S turns out to be 

S' = -q- S 2 . 

We assume that the differential equation contains a parameter v and that we are inter- 
ested in the large-z/ asymptotic of S. The particular choice of the asymptotic expansions 
below is the relevant case for the consideration of the asymptotic of Bessel functions. 
We assume that as v — > oo, the function q has the asymptotic expansion 
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In that case, the function S can be seen to have the asymptotic form 

oo 



i=-l 



Using these asymptotic forms in the differential equation for S, the asymptotic orders 
are seen to be given by 



S-i = ±V-<?-2, 

Let us apply these formulas to the differential equation for the Bessel functions I u+a (i/z) 
and K v+a (i/z). The relevant differential equation reads 

Therefore, for the given example, we find 

Linearly independent solutions of the differential equation are proportional to I u+a (uz) 
and K v+a (vz). From the known behaviour of these functions for large arguments [26J, 
we can conclude that >S_i = q~2 corresponds to I u + a , whereas 5_i = —y/—q-2 

corresponds to K v+a . For the present occasion we need the asymptotics for I v+a and 
continue with this case only. Using the asymptotic orders in S, the solution has the 
asymptotic behavior 



z) ~ const exp < v \ dzS^i > exp < dz ( S { 



?0 " i)} exp {X^ */ dzS ^j 



The constant prefactor is determined from the known Debye expansion of the Bessel 
function I v {yz) and it reads 

1 

const = . 

Using the explicit form of g_2 and g_i for the present example, one obtains therefore the 
result 

a/2 



A/27 
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where t = 1/VTTz 2 , (3 = VTTz 2 + ln[z/(l + y/TTz 2 )}, 



Ci(z,a) = J dzS h 

and with the understanding that for a = the known answers for l v {yz) are reproduced 
[35] . The list of the first Ci(z,a) obtained are 

t a 2 t at 2 5t 3 
ClM = 8~ — ~ — "24' 

„ , . t 2 a 2 t 2 13at 3 a 3 t 3 3t 4 a 2 t 4 5at 5 5t 6 

C 2 (z,a) = 1 1 1 1 , 

2V ' ; 16 4 24 6 8 2 8 16 ' 

25t 3 13a 2 t 3 aH 3 7at 4 a 3 t 4 531t 5 7a 2 t b aH 5 

C 3 (z,a) = 1 1 1 

iK ' ; 384 48 24 8 2 640 4 8 

Hat 6 2aH 6 221t 7 25a 2 t 7 15at 8 1105t 9 . . 

+ 1 . (45) 

4 3 128 16 8 1152 v ; 

Higher orders can be produced as needed. 
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